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ABSTRACT 

Exoplanetary systems closest to the Sun, with the brightest host stars, provide the most favorable 
opportunities for characterization studies of the host star and their planet (s). The Transit Ephemeris 
Refinement and Monitoring Survey uses both new radial velocity measurements and photometry in 
order to greatly improve planetary orbit uncertainties and the fundamental properties of the star, 
in this case HD 130322. The only companion, HD 130322b, orbits in a relatively circular orbit, 
e = 0.029 every ^10.7 days. Radial velocity measurements from multiple sources, including 12 
unpublished from the Keck I telescope, over the course of ^14 years have reduced the uncertainty 
in the transit midpoint to ^2 hours. The transit probability for the 6-companion is 4.7%, where 
MpSini = 1.15Mj and a = 0.0925 AU. In this paper, we compile photometric data from the TII 
0.8m Automated Photoelectric Telescope at Fairborn Observatory taken over ^14 years, including 
the constrained transit window, which results in a dispositive null result for both full transit exclusion 
of HD 130322b to a depth of 0.017 mag and grazing transit exclusion to a depth of ^O.OOI mag. Our 
analysis of the starspot activity via the photometric data reveals a highly accurate stellar rotation 
period: 26.53±0.70 days. In addition, the brightness of the host with respect to the comparison stars 
is anti-correlated with the Ca H H and K indices, typical for a young solar-type star. 

Subject headings: planetary systems - techniques: photometric - techniques: radial velocities - stars: 
individual (HD 130322) 


1. INTRODUCTION 

Studying the ephemerides, or orbital parameters, of 
nearby planets is one of the oldest sub-helds in as¬ 
tronomy. It took nearly 1500 years for a new celestial 
model to supplant Ptolemy’s stationary, geocentric sys¬ 
tem specifically because, through the use of a number 
of clever maneuvers (equants, epicycles, and deferents), 
he was able to accu rately predict th e motion of the So¬ 
lar System planets (|Gingerichl II997I1 . And it was pre¬ 
cisely because Copernicus’ heliocentric model did not 
predict accurate planetary phenomena that nearly 100 
years passed before the theory was generally accepted 
by the scientihc community. With the purpose of im¬ 
proving orbital uncertainties and fundamental properties 
of the host star, the Transit Ephemeris Refinement and 
Monitoring Survey (TERMS) team s eeks to characteriz e 
individual nearby planetary systems (jKane et al.l 12009(1 . 
The giant planet around HD 130322, a KO star, was 
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first detected by iTJdrv et ^ (120001 1 using the radial ve¬ 
locity (RV) technique via the CORALIE echelle spectro¬ 
graph. They reported the period of the planet was P = 
10.720 ± 0.007 days, m sinj = 1.02 Mj, and eccentricity 
of e = 0.044 ± 0.01 8. The planet was later confirmed by 
iButler et al.l (|2006ll . who observed an additional 12 RV 
measurements using Keck and determined that the pe¬ 
riod was P = 10.70875 ± 0.00094 days, m sini = 1.089 
Mj, and an eccentricity of e = 0.025 ± 0.032. Th e giant 
planet was further observed bv lWittenmver et al.l (|2009f ) 
using both the Hobby-Eberly Telescope (HET) and the 
Harlan J. Smith telescope. They found the combina¬ 
tion of all four datasets produced a large root-mean- 
square (rms) variability of 14.8 ms“*^ and an anoma¬ 
lous periodicity at 35 days, specifically due to CORALIE 
data. Th eir orbital parame ters, without the RV measure¬ 
ments bv lUdrv et al.l (|2000ll . found P — 10.7085±0.0003, 
m sini = 1.04 ± 0.03 Mj, and e = O-O l l 0 .020. Be¬ 
cause of the small eccentricity, iTrillingl (|2000[ 1 was not 
able to put a lower bound on the mass of the planet as 
determined by tidal constraints, only an upper limit of 
43.8 Mj. O bservations using the Sp i tzer in frared spec¬ 
trograph by iDodson-Robinson et al.l (|2011h resulted in 
the detection of a debris disk around the host star. 

Here we present our complete RV data set from a 
number of sources (including those mentioned above in 
addition to previously unpublished measurements from 
the 10m Keck I telescope) that has a time baseline of 
^^14 years, discussed in (JS) The analysis of the Keple- 
rian orbital solution of these data produce refined orbital 
ephemeris for the host star HD 130322, with a predicted 
transit depth of 1.57%, and Icr transit window of 0.329 
days (lj3|). In SjU we determine the differential magni¬ 
tude of the host star with respect to multiple comparison 
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stars in order to better understand seasonal and nightly 
brightness fluctuations. The evaluation of starspot vari¬ 
ability f H4.1|) allows us to calculate a stellar rotation pe¬ 
riod, while understanding the stellar magnetic activity 
( M.2i gives us insight into the age of HD 130322. 

2. HOST STAR PROPERTIES 

The HD 130322 system has been monitored via the 
RV technique several times in the past. We provide 
a complete RV data set that consists of 118 measure¬ 
ments acquired with CORALIE at the 1.2m Euler- 
Swiss telescope (|Udrv et al.ll2000ll . 35 measurements ac¬ 
quired with the 2.7m Harlan J. Smith telescope and 
the High Resolution Sp ectrograph (HRS) on the HET 
(jWittenmver et all2009H , and 24 measurements acquired 
with the HI RES echelle spectr ometer on the 10m Keck 
1 telescope (|Butler et al.M2006ll . the most recent 12 of 
which are previously unpublishedl3 We use this com¬ 
bined data set (shown in Table [IJ in order to calculate 
the fundamental properties of the star as well as the Ke- 
plerian orbital solution of the planet. 

2.1. Fundamental Parameters 

We derive the host star properties by fitting the 
high resolution HIRES, HR S, and CORALIE da t a wit h 
Spectroscopy Made Easy (jValenti fc Piskuno^ 1199611 . 
or SME, via the wavelength in t ervals , line data, and 
technique of I Valenti fc Fischeil (l2005l). We app lied 
the revised Hipya r cos pa rallaxes (jvan Leenwenl 12?)?)^ to 
the I Valenti et al.l (|20n9[ l methodology as well as sur- 
face gravity from a Yq nsei-Yale stellar structure model 
((Demaraue et al.l 1200411 . As a result, the fundamen¬ 
tal stellar parameters are listed in Table [2 where V- 
magnitude and distance were determined by Hipparcos, 
B — V from Tycho-2, while effective temperature, sur¬ 
face gravity, projected rotational velocity, stellar mass, 
and stellar radius were a result of SME. The high pre¬ 
cision of the stellar radius, namely R* = 0.85±0.04 Rq, 
is important when determining the depth and duration 
of a possible planetary transit. As a comparison to the 
SME result, the stellar radius was determined using the 
Torres relation ((Torres et al.ll20ldll : R* = 0.88±0.04 Rq. 
We have also conducted an empirica l surface brightness 
calculation oer lBovaiian et all (1201411 which averages the 
V — J, V — P[, and V — K surface brightness relations, 
resulting in an angular diameter of 0.252 ± 0.006 mas. 
Folding in the parallax and error, the radius is 0.89±0.04 
Rq. All three of these techniques show a very strong 
consensus for both the stellar radius and error of the 
host star. The iron abundance, [Fe/H], as determined 
by SME, as well as other element abundances, will be 
discussed in (12.21 

W hile our results are co nsistent with a typical K-type 
star (|Bovaiian et al.ll2ni2[ l. we note th e diffe rences be¬ 
tween stellar properties in lUdrv etldl (|2000ri . namely, 
their Table 1 an d our Table [2 We have used the up¬ 
dated Hipparcos (Ivan Leeuwenll200^ catalog, which may 
account for the varying B — V and distance determina¬ 
tions; they cited B — V = 0.781 and the distance to be 

W Based on observations obtained at the W. M. Keck Observa- 
tory, which is operated jointly by the University of California and 
the California Institute of Technology. Keck time has been granted 
by the University of Hawaii, the University of California, Caltech, 
and NASA. 


29.76 pc. Our effective temperature, 5387 ± 44, is also 
-|-50 K above their referenced temperatures. 

2.2. Stellar Abundances 

Stellar abundances have been measured for HD 130322 
by at least a dozen differ en t groups, f or ex ¬ 
ample iV alenti fc Fischeil (1200511 : iNeves et al.l (|2009ll : 
iDelsrado Mena et al.l ( 2(Tlfll l. Due to the proximity of 
the host star to the Sun, 31.54 pc, a wide variety of 
elements have been measured within HD 130322, from 
a-type to neutron-captur e. Using the same analysis as 
seen in Hypatia Catalog (|Hinkel et al.ll20i2l . we renor¬ 
malized the abundanc e measurements for e ach dataset to 
the same solar scale (jLodders et al.l[2(i0^ and then de¬ 
termined the maximum measurement variation between 
the groups, or the spread^ to quantify the consistency 
of the abundances. When analyzing the abundances in 
the Hypatia Catalog, element measurements were only 
considered where the spread was less than the respec¬ 
tive error bar associated with that element, or where 
group-to-group variations were small, and then the me¬ 
dian value was used. For HD 130322, [Fe/H] = 0.12 dex, 
however the spread between the groups was 0.23 dex. In 
other words, the iron ratio was not agr eed upon by the 
variou s groups, where the renormalized lEciivillon et al.l 
(l2Q114f) measurement w as [Fe/H] = 0.04 and renormalized 
iBodaghee et al.l (j2003f l determined [Fe/H] = 0.27 dex. 
For this reason, HD 130322 was not included in the anal¬ 
ysis (or reduced version) of the Hypatia Catalog. Per the 
SME analysis, [Fe/H] = 0.07 ± 0.03 dex, using the dis¬ 
cussio n and solar abundance scale in IValenti fc Fischeil 
(l2005fl . The renormalization gives [Fe/H] = 0.11, which 
is very close to the median value found for the other data 
sets in Hypatia. 

There were a number of other elements within the star 
that were measured by different groups. Per the Hypatia 
analysis, where the spread in the abundances were less 
than respective error and the median value taken: [N/Fe] 
= -0.14 dex, [Al/Fe] = -0.12 dex, [S/Fe] = -0.14 dex, 
[Cu/Fe] = -0.15 dex, [Sr/Fe] = 0.07 dex, [YH/Fe] = -0.09 
dex, [Ball/Fe] = -0.06 dex, Ce/Fe] = 0.07 dex, [Cell/Fe] 
= -0.05 dex, and [EuH/Fe] = -0.18 dex. In general, we 
find that the majority of abundances well measured in 
HD 130322 are significantly sub-solar. Despite measuring 
[Fe/H] = 0.07 ± 0.03 via SME, not much can be said 
conclusively about the overall [Fe/H] content, given the 
large spread in the abundances determined by different 
methods. 

3. KEPLERIAN ORBIT AND TRANSIT EPHEMERIS 

We fit a Keplerian orbital solution to the RV data 
(shown in Table [T]) using the partiall y linearized, least- 
square s fitting procedure described in lWright fc Howard 
(j2009fl with parameter uncertainties estimat ed using the 
BOOT TRAN bootstrapping routines from iWang et al.l 
(|2012f) . The resulting Keplerian orbital solution is shown 
in Table [3] where the stellar parameters for the host star 
described in Section [2] and summarized in Table [2] were 
used to determine the minimum mass and semi-major 
axis of the planet. The phased data and residuals to the 
fit are shown in Fig. [TJ We find the offsets to be 24.3, 
24.7, -27.2, and -23.6 ms“^ for data from 2.7m McDonald 
Observatory telescope (2.7m), CORALIE, HIRES (pre- 
2004 or BJD prior to 13005.5), and HIRES (post-2004), 
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respectively. Regarding the CORALIE data, the me¬ 
dian velocity value was subtracted from the data and 
the velocities were converted from km/s to m/s, which 
were then used to calculate the offsets with respect to 
the HET’s HRS data. The fit including the CORALIE 
data has a x^ed = 1-35 and RMS = 14.6 ms“^. With¬ 
out the CORALIE data, these numbers change to 1.46 
and 8.67, respectively. However, the time baseline of the 
CORALIE data signihcantly improves the determination 
of the orbital period so our fit includes these data for the 
subsequent analysis. 


TABLE 1 

Radial velocities measured eor HD 130322 


BJD 

(-2440000) 

RV 

(m/s) 

zblfT 

(m/s) 

Tel 

11755.76855 

-54.8 

1.0 

HIRES 

11984.06010 

-102.7 

1.4 

HIRES 

12065.94426 

-55.6 

1.4 

HIRES 

12127.80910 

77.8 

1.3 

HIRES 

12128.76352 

44.8 

1.4 

HIRES 

12162.72653 

-60.7 

1.4 

HIRES 

12335.12452 

-127.8 

1.3 

HIRES 

12488.77933 

-24.9 

1.5 

HIRES 

12683.09362 

54.9 

1.4 

HIRES 

12805.91934 

-94.5 

1.4 

HIRES 

13153.85287 

42.5 

1.3 

HIRES 

13426.11560 

-43.4 

1.1 

HIRES 

13842.00573 

48.6 

1.1 

HIRES 

15351.81672 

75.5 

1.1 

HIRES 

15636.05652 

-90.5 

1.1 

HIRES 

15673.83685 

25.2 

1.2 

HIRES 

15700.79799 

-59.4 

1.5 

HIRES 

15700.80123 

-62.7 

1.3 

HIRES 

15734.89254 

56.2 

1.2 

HIRES 

15789.75142 

79.1 

1.3 

HIRES 

15961.16166 

86.0 

1.2 

HIRES 

16000.02899 

-102.0 

1.2 

HIRES 

16075.79698 

-51.6 

1.2 

HIRES 

16451.82279 

23.7 

1.2 

HIRES 

13585.64900 

83.7 

7.5 

2.7m 

13843.89253 

-18.0 

7.5 

2.7m 

13863.78301 

75.5 

8.6 

2.7m 

13910.78043 

-68.5 

8.1 

2.7m 

14251.84318 

-72.8 

9.4 

2.7m 

13471.80558 

-99.8 

7.2 

HRS 

13481.88526 

-106.9 

6.6 

HRS 

13486.85864 

105.1 

6.4 

HRS 

13488.75815 

72.2 

5.9 

HRS 

13509.79117 

101.3 

6.1 

HRS 

13512.78123 

-65.6 

5.1 

HRS 

13527.74971 

27.0 

6.2 

HRS 

13542.69985 

55.4 

5.6 

HRS 

13543.70614 

-4.9 

6.1 

HRS 

13550.70420 

105.5 

6.1 

HRS 

13837.89677 

-12.6 

5.9 

HRS 

13842.88880 

27.6 

6.3 

HRS 

13868.80896 

-78.8 

5.7 

HRS 

13882.78043 

83.0 

6.0 

HRS 

13897.72683 

-44.2 

6.1 

HRS 

13900.72079 

-85.4 

6.0 

HRS 

13936.63557 

110.4 

6.6 

HRS 

14122.01834 

-12.7 

6.8 

HRS 

14128.00335 

47.0 

6.7 

HRS 

14135.98084 

-113.2 

6.5 

HRS 

14139.97029 

89.9 

7.2 

HRS 

14140.96840 

98.1 

6.1 

HRS 

14144.96962 

-99.4 

6.6 

HRS 

14157.01611 

-112.4 

6.8 

HRS 

14158.92425 

-40.9 

6.8 

HRS 

14163.92465 

26.3 

6.7 

HRS 

14168.90656 

-71.6 

6.5 

HRS 

14173.98269 

69.7 

7.4 

HRS 

14176.87914 

-90.6 

5.7 

HRS 


TABLE 1 — Continued 


BJD 

RV 

dzlfj 

Tel 

(-2440000) 

(m/s) 

(m/s) 



14191.92631 

11257.85195 

11267.80486 

11267.81699 

11273.86307 

11287.71123 

11287.72337 

11291.77019 

11294.82378 

11295.76291 

11296.60937 

11296.83574 

11297.61667 

11297.82760 

11298.60898 

11298.83006 

11299.60924 

11299.82827 

11300.60428 

11300.82617 

11301.59933 

11301.81846 

11302.73856 

11304.74462 

11305.79957 

11306.70099 

11307.78007 

11307.79214 

11308.77008 

11308.78216 

11309.74470 

11309.75677 

11310.65574 

11310.66781 

11311.76937 

11312.71231 

11313.70809 

11314.71380 

11315.65132 

11316.70334 

11317.73465 

11318.71336 

11319.67233 

11320.73422 

11320.74629 

11335.68044 

11335.69251 

11336.59905 

11336.61110 

11339.65310 

11339.66512 

11340.65091 

11340.66296 

11342.62586 

11342.63788 

11355.64903 

11364.57708 

11366.58262 

11367.59944 

11368.57076 

11369.60470 

11369.61675 

11370.51479 

11370.52684 

11373.52546 

11373.53762 

11374.54843 

11374.56049 

11375.55068 

11375.56276 

11376.51940 

11376.53149 

11380.59301 

11380.60510 

11381.51359 

11381.52568 


20.5 

6.2 

99.0 

9.0 

39.0 

9.0 

45.0 

9.0 

2.0 

9.0 

-52.0 

9.0 

-68.0 

9.0 

137.0 

9.0 

44 

15 

-31.0 

9.0 

-72 

10 

-86 

10 

-86 

10 

-99 

10 

-71 

10 

-52 

10 

-18.0 

9.0 

2 

10 

53.0 

9.0 

69 

10 

103 

10 

119.0 

9.0 

105 

12 

71 

11 

-1 

10 

-58 

10 

-129 

14 

-122 

14 

-131 

10 

-134 

10 

-52 

11 

-47 

10 

13 

10 

5 

10 

51 

34 

130.0 

9.0 

123.0 

9.0 

100 

12 

46 

24 

-21 

14 

-81 

11 

-95 

13 

-73.0 

9.0 

-37 

11 

-40 

11 

141 

14 

122 

14 

73 

14 

84 

13 

-100 

17 

-122 

17 

-109 

12 

-109 

12 

-16 

11 

0 

10 

119 

10 

43 

11 

146 

11 

124 

11 

83 

11 

41 

15 

49 

15 

-7 

15 

-46 

17 

-66 

27 

-103 

27 

2 

18 

-14 

17 

47 

11 

57 

11 

58 

21 

71 

19 

19 

15 

-14 

16 

-46.0 

9.0 

-47.0 

9.0 


HRS 

CORALIE 

CORALIE 

CORALIE 

CORALIE 

CORALIE 

CORALIE 

CORALIE 

CORALIE 

CORALIE 

CORALIE 

CORALIE 

CORALIE 

CORALIE 

CORALIE 

CORALIE 

CORALIE 

CORALIE 

CORALIE 

CORALIE 

CORALIE 

CORALIE 

CORALIE 

CORALIE 

CORALIE 

CORALIE 

CORALIE 

CORALIE 

CORALIE 

CORALIE 

CORALIE 

CORALIE 

CORALIE 

CORALIE 

CORALIE 

CORALIE 

CORALIE 

CORALIE 

CORALIE 

CORALIE 

CORALIE 

CORALIE 

CORALIE 

CORALIE 

CORALIE 

CORALIE 

CORALIE 

CORALIE 

CORALIE 

CORALIE 

CORALIE 

CORALIE 

CORALIE 

CORALIE 

CORALIE 

CORALIE 

CORALIE 

CORALIE 

CORALIE 

CORALIE 

CORALIE 

CORALIE 

CORALIE 

CORALIE 

CORALIE 

CORALIE 

CORALIE 

CORALIE 

CORALIE 

CORALIE 

CORALIE 

CORALIE 

CORALIE 

CORALIE 

CORALIE 

CORALIE 
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TABLE 1 — Continued 


TABLE 3 — Continued 


BJD 

(-2440000) 

RV 

(m/s) 

zblfT 

(m/s) 

Tel 

11382.48455 

-85 

10 

CORALIE 

11382.49669 

-78 

10 

CORALIE 

11382.56236 

-83 

11 

CORALIE 

11382.57447 

-87 

11 

CORALIE 

11383.53612 

-71 

10 

CORALIE 

11383.54681 

-88 

10 

CORALIE 

11384.54932 

-71.0 

9.0 

CORALIE 

11384.56142 

-62 

10 

CORALIE 

11385.54622 

12 

10 

CORALIE 

11385.55834 

10 

11 

CORALIE 

11386.54337 

86 

10 

CORALIE 

11386.55545 

74 

11 

CORALIE 

11388.53150 

132 

14 

CORALIE 

11388.54362 

145 

14 

CORALIE 

11389.49681 

127 

14 

CORALIE 

11389.50891 

121 

12 

CORALIE 

11390.46734 

73 

11 

CORALIE 

11390.47944 

81 

11 

CORALIE 

11391.46814 

5.0 

9.0 

CORALIE 

11391.48030 

11.0 

9.0 

CORALIE 

11392.52992 

-48.0 

9.0 

CORALIE 

11392.54196 

-38 

10 

CORALIE 

11393.51766 

-80 

10 

CORALIE 

11393.52971 

-91 

10 

CORALIE 

11394.47208 

-68 

11 

CORALIE 

11394.48410 

-78 

11 

CORALIE 

11395.47053 

-25 

17 

CORALIE 

11395.48250 

-42 

11 

CORALIE 

11397.47145 

88 

20 

CORALIE 

11397.48921 

90 

15 

CORALIE 

11398.47062 

128.0 

9.0 

CORALIE 

11398.48258 

121.0 

9.0 

CORALIE 

11399.47201 

139 

10 

CORALIE 

11400.48982 

118 

10 

CORALIE 

11401.46724 

63 

11 

CORALIE 

11402.48058 

-1 

17 

CORALIE 

11403.47387 

-52 

12 

CORALIE 

11404.47239 

-109 

12 

CORALIE 

11405.47319 

-46 

22 

CORALIE 

11406.48267 

-6 

21 

CORALIE 

11412.48141 

28 

11 

CORALIE 

11412.49345 

35 

11 

CORALIE 

11424.49242 

-61 

20 

CORALIE 


Note. — The line separates the pre-2004 and post-2004 HIRES 
data. 


TABLE 2 

Stellar Parameters 


Parameter 

Value 

Source 

V 

8.04 

Hipparcos 

B-V 

-0.16 

Tycho-2 

Distance (pc) 

31.54 ±1.18 

Hipparcos 

na (K) 

5387 ± 44 

SME 

logg 

4.52 ±0.06 

SME 

t/sinz (kms“^) 

0.5 ±0.5 

SME 

M, (Mq) 

0.92 ±0.03 

SME 

R* {Rq) 

0.85 ±0.04 

SME 


TABLE 3 

Keplerian Fit Parameters 


Parameter Value 


P (days) 10.70871 ± 0.00018 

Tc “ (JD - 2440000) 16745.594 ± 0.085 

Tp (JD - 2440000) 13996.4 ±1.1 

e 0.029 ±0.016 


Parameter 

Value 

K (ms“l) 

112.5 ± 2.4 

LO (deg) 

193 ± 36 

X?ed 

1.35 

RMS (ms-l) 

14.60 


^ Time of transit. 

^ Time of periastron passage. 

The lack of a linear trend over a long period time in 
Fig. [1] opens up the possibility to constrain the presence 
of additional companions in the HD 130322 system. If 
msin* is the “minimum mass” (accounting for inclina¬ 
tion) then we can also consider the minimum value that 
m sin i could possibly have given a linear trend that per¬ 
sists over time (the “ minimum minimu m mass”), which 
is described in detail Ihhng et al.l (I2015D . Given that we 
only have an upper limit on a trend, we have measured 
the maximum value that the minimum msini could take, 
or the “maximum minimum minimum mass” {Mmmm). 
We use BOOTTRAN to find the Ict maximal value of the 
linear velocity over years, where \dv/dt\ = 0.0047 
m/s/day = 1.716 m/s/yr. In additi on to the values fr om 
Tabled! we employ Equation 1 in iFeng et al.l (I2015D to 
calculate Mmmm = 1.83Mj as an upper-bound for a 
possible additional companion in the HD 130322 system. 

Finally, we use the revised orbital properties of the 
planet described above to derive the predicted transit 
properties. The predicted time of mid-transit produced 
by the Keplerian orbital solution is Tc=16745.594±0.085 
(see Table [3|). Since the orbit is close to circular, the ec- 
centricity has a negligible effect on the transit properties 
(|Kane fc von Bra unll2008D. Using the mass-radius rela- 
tionship of iKane fc Gelinol (120121) . we adopt a radius of 
the planet of Rp = l.ORj. These combined parameters 
for the planet result in a transit probability of 4.7%, a 
predicted transit duration of 0.16 days, and a predicted 
transit depth of of 1.57%. The size of the Icr transit win¬ 
dow is 0.329 days which can b e adequately mon itored in 
a single night of observations (iKane et al.ll2009D . 

4. PHOTOMETRY 

Between Jan 2, 2001 and June 26, 2014, 1569 observa¬ 
tions were obtained for HD 130322 at Fairborn Observa¬ 
tory in Arizona using the Til 0.80m APT. The APT is 
able to determine the differential brightness of the pri¬ 
mary star HD 130322 (P: V = 8.04, B-V = 0.781, KOV) 
with respect to three comparison stars: HD 130557 (Gl: 
V=6.15, B-V=-0.02, B9V), HD 129755 (G2: V=7.58, B- 
V=0.41, F2), and HD 132932 (G3: V=7.74, B-V=0.40, 
F2). In the initial 2001 observing season, we found that 
our original comparison star C3 was a low-amplitude 
variable, so we replaced it with HD 132932 in 2002. 
Therefore, we have only two comparison stars, Cl and 
C2, in common for all 14 observing seasons, whereas sea¬ 
sons 2-14 have in common the three comparison stars 
given above. Like the other telescopes operated on site 
by Tennessee State University, the Stromgren b and y 
bands are separated and concurrently measured by a 
photometer with two-channel precision, a dichroic fil- 
ter, and two EMI 9124QB bi-alkali photomultiplier tubes 
(|Henrvlll999f) . 

We compute the six permutations of the differential 
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Fig. 1. — The Keplerian orbital solution using all of the data shown in Table [T] resulting in the fit parameters shown in Table [S] The 
typical internal error bars for each data point are plotted. Left: RV data phased on the best-fit solution, where the origin of the data is 
indicated by the different symbols, shown in the figure. Right: Residual velocities with respect to the fitted orbital solution, with the same 
symbols as the left-side panel. 


magnitudes of the four stars in a combinitorical fashion, 
namely P - Cl, P - C2, P - C3, (73 - Cl, (73 - (72, 
and (72 — (71. The magnitudes are then corrected for ex¬ 
tinction due to the atmosphere and transformed to the 
Stromgren system, such that the differential b and y ob¬ 
servations are combined into a single {b + y)/2 band, 
indicated with the subscript by. To achieve the maxi¬ 
mum possible precision, we also combine the three com¬ 
parison stars to determine the differential magnitudes of 
HD 130322 with respect to the mean brightness of the 
comparison stars. The precision of the individual dif¬ 
ferential magnitudes P — ((71 + C2 + (73)/3by ranges be¬ 
tween 0.0010 mag and ^ 0.0015 mag on clearer nights, 
as determined from the nightly scatter o f the compar i¬ 
son stars. Further details can be found in iHenrvI (jl999H . 
lEaton et al.l (|2003l) . and references therein. 

The 1470 individual P — {Cl PC2P Ci)liby differen¬ 
tial magnitudes computed from the 13 observing seasons 
(2002-2014) are plotted in the top panel of Fig. [2l The 
observations are normalized so that all 13 seasons have 
the same mean as the hrst season 2002, indicated by 
the horizontal line in the top panel, to remove season- 
to-season variability in HD 130322 caused by a possible 
starspot cycle (see below). The normalized nightly ob¬ 
servations scatter about their grand mean of 1.02748 mag 
with a standard deviation of ct = 0.00331 mag, which 
is more than a factor of two larger than the ^ 0.0010 
- ~ 0.0015 mag measurement precision, which suggests 
HD 130322 has nightly low-amplitude variation. 

The normalized differential magnitudes from the last 
13 observing seasons are shown in Fig. [2] (middle panel), 
where they are phased with the planetary 10.7 day or¬ 
bital period and the mid transit time (Tc) given in Table 
[21 A fit using a least-squares sinusoid provides a photo¬ 
metric semi-amplitude of 0.00023 ± 0.00011 mag and 
places a one milli-magnitude (0.001 mag) upper bound 
on the brightness variability of th e host star. In add i¬ 
tion, per simi l ar res ults fo und in iQue oz et al.l (j2001l ). 
iPaulson et al.l (j2(104[l . and iBoisse et al. (j2012l )7 we dis¬ 
miss the possibility that jitter induced stellar activity 
can account for the 10.7-day RV fluctuations. The con¬ 


stancy of the photometric measurements reveals the true 
planetary reflex motion seen in the RV variations of HD 
130322 is a result of the orbiting planet. 

A closer view of the predicted transit window is shown 
in the bottom panel of Fig. |2l plotted with an expanded 
abscissa. Similar to the middle panel, the solid curves 
shows the predicted central transit, phased at 0.0, for 
a duration of 0.16 days or ~ 0.015 units of phase and 
a depth of 1.57% or ~ 0.017 mag. These values were 
determined using the stellar radius (Table [5]) and or¬ 
bital ephemeris of the planet (Table [3]). The ±lcr un¬ 
certainty in the transit window timing, as determined 
by the error bars for both the stellar radius (Table [2]) 
and the improved orbital ephemeris (Table jS]), is in¬ 
dicated by the vertical dotted lines. There are 1405 
observations that lie outside the predicted transit win¬ 
dow, which have a mean of 1.027490 ± 0.000089 mag. 
The 65 observations that fell within the transit window 
have a mean of 1.027278 ± 0.000303 mag. The dif¬ 
ference in these two light levels is our “observed transit 
depth,” —0.00021 ± 0.00032 mag, consistent to four 
decimal places. Thus, we are able to rule out full tran¬ 
sits, a dispositive null result, with a predicted depth near 
0.017 mag and also grazing transits near the predicted 
time to a depth of ^0.001 mag. 

4.1. Starspot Analysis 

Given that the scatter of 0.00331 mag in the normal¬ 
ized data set of Fig. [1] is significantly larger than the 
observation precision, we suspect low-amplitude, night- 
to-night starspot variability in HD 130322. Inspection 
of the top panel of Fig. [5] reveals differences in the 
amount of scatter from year to year that could also be 
caused by starspot activity. A solar-type star’s rotation 
period may be determined from the rotational modula¬ 
tion of starspots on the photosp here by measur i ng th e 
variation in stellar brightness per iSimoson et al.l ()20inH . 
In addition, starspots can resemble an or biting planet 
by g enerating periodic RV fluctuations (iQueloz et al.l 
I 2 OOII) . Therefore, to determine the behavior of poten¬ 
tial starspots, we analyzed all 13 seasons of normalized 
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Fig. 2.— Top: the 1470 differential magnitudes (P - (Cl + C2 + 
C3) / ^by) for HD 130322, taken using the 0.8m APT from 2002- 
2014, where all 13 observing seasons are normalized to standardize 
the yearly average. Middle: Observations phased to the planet’s 
ephemeris. The orbital phase curve semi-amplitude is 0.00023 zb 
0.00011 mag, fit with a least-squares sine, which confirms the b- 
planet with the lack of periodic light variability in the host star. 
Bottom: A zoomed-in portion of the middle-plot, centered on the 
central transit midpoint. The solid curve shows the predicted cen¬ 
tral transit, with a depth of 1.57% or 0.017 mags and duration of 
0.16 days or 0.015 units of phase. The vertical dotted lines show 
the dilcr transit window. Transits are excluded to a depth of 0.017 
mag while grazing transits are excluded to a depth of ^^O.OOl mag. 


Trial Frequency (c/d) 



Photometric Phase 

Fig. 3.— Top: the 91 differential magnitudes (P - (Cl -h C2 -b 
C3) / ^by) for HD 130322, taken during the 2013 season. Middle: 
A frequency spectrum of the 2013 observing season of HD 130322, 
such that the best frequency is at 0.03773 ib 0.00036 cycles per day 
(or c/d). Bottom: the 2013 seasonal observations phased with the 
corresponding best rotational period of 26.53 ib 0.70 days. The 
peak-to-peak amplitude of 0.011 mag shows coherent variability, 
which may be due to the rotational modulation of starspots, and 
low noise around the trend. All 14 observing seasons exhibit similar 
modulation (see Table 0. 


photometry using a periodogram analysis. While low- 
amplitude (0.002-0.017 mag) periodic brightness fluctu¬ 
ations were found during each season, there was no un¬ 
usual periodicity for the yearly C2 —Cl comparison stars. 
The frequency spectrum for the penultimate 2013 observ¬ 
ing season is shown in the top panel of Figure [U while 
the phase curve is given in the bottom panel. 

Results of our complete seasonal period analyses are 
given in Table 21 Here we include the first observing 
season from 2001, in which comparison star 3 was later 
found to be variable and was replaced for the subsequent 
seasons, as previously mentioned. The period analysis of 
season 1 (2001) is therefore based on differential magni¬ 
tudes computed as P — (Cl -|- C2)l2by. Another minor 
caveat should be noted for the 2004 and 2011 observ¬ 
ing seasons. In both cases, the first harmonic of the 
rotation period was found to have a slightly higher peak 
than the rotation period. For these two seasons, we es¬ 
timated the rotation periods and their uncertainties by 
doubling both values. The mean of the 12 rotation pe¬ 
riods, excluding 2004 and 2011, is 26.53 days. The in¬ 
dividual rotation periods scatter about that mean with 
a standard deviation of 2.44 days, and the standard de¬ 


viation of the mean is 0.70 days. Therefore, we deter¬ 
mine that 26.53 ± 0.70 days is our most accurate cal¬ 
culation for the star’s rotation period, which matches 
well w ith the rotation period derived by iSimpson et al.l 
(IMoli . or 26.1 ± 3.5 days, from the first six seasons of our 
Apt data set. We note the significantly smaller uncer¬ 
tainty in our determination due to the greatly extended 
14 year baseline. We found typical starspot-filling fac¬ 
tors of one percent or less in Column 6 of Table 01 where 
the peak-to-peak amplitudes for each season range from 
0.002-0.017 mag. Because the stellar rotation period of 
26.53 days is more than a factor of two from the 10.7 day 
planetary RV period and harmonics, in addition to the 
small impact from starspots, we find that the 10.7-day 
planet could not be the result of stellar activity. 

Examining the stellar rotation with respect to incli¬ 
nation, we first assume that the inclination of rotation 
axis for HD 130322 is close to 90°, in which case the 
stellar radius (0.85 Mq) and nsinz (0.5 km s“^) predict 
a stellar rotation period of ~ 85.5 days. Because this 
value is over three times longer than our observed value 
of Prot = 26.53 days, the logical conclusion is that the 
stellar rotation axis must have a low inclination and a 
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low planetary transit probability. However, if we substi¬ 
tute t he value of fsini = 1.61 km s“^ from iButler et all 
(j2006ll . the predicted rotation period is 26.55 days, iden¬ 
tical to our observed rotation period within the uncer¬ 
tainties. This implies a very high inclination near 90° 
and, therefore, a high probability of transits. Unfortu¬ 
nately, our photometry unambiguously rules them out. 

4.2. Magnetic Activity 

To look for evidence of magnetic cycles in HD 130322, 
we analyze the variability in the Ca H H and K indices, 
both proxies for stellar magn etic activity (|Baliunas et al.l 
119951: [Lockwood et al.ll200^ . and APT photometry over 
the entire observing period. These magnetic cycles could 
potentially resemble the period of a long-period planet 
within the RV data. At the top of Fig. 21 we show 
the seasonal means for the Mount Wilson S-index as 
determined f r om th e Keck I RV spectr a , desc ribed in 
IWright et al.l (|2004[ 1 ; [Isaacson fc Fischt^ (j201fl[ l . While 
we do not have the Keck H and K (or RV) measure¬ 
ments for all 13 of our photometric observing seasons, 
there is notable variability on the order of several years. 
The middle two panels show the seasonal mean of HD 
130322 (P) varying with respect to the two comparison 
stars (Cl and C2) throughout the 14-year observations 
without normalization (Table 21) ■ The horizontal dot¬ 
ted line indicates the standard deviation of each seasonal 
mean as compared to the grand mean, given numerically 
in the lower right corner. The range in magnitude of the 
seasonal mean is printed in the lower left corner. The 
brightness curves in the middle of Fig. Ulshow HD 130322 
varying on the order of multiple mmag with respect to 
both of the comparison stars. Similarly, the yearly av¬ 
erage of the comparison stars (72 — (71 is given in the 
bottom panel, with a very small standard deviation of 
0.0005 mag. Since the comparison stars demonstrate 
photometric stability over the 13 observing seasons, the 
fluctuations seen in middle two panels must be intrinsic 
to the host-star HD 130322. 

The variations in H and K and APT observations plot¬ 
ted in Fig. 21 appear to be cyclic. Analyses of the yearly 
means for the Ca H indices, P — (71, and P — (72 using 
a least-squares, sine-fit periodogram results in the same 
periods to within uncertainty: 5.22 ± 0.16, 5.19 ± 0.20, 
and 5.1208 ± 0.22 yr, respectively. These amplitudes 
and timescales for HD 130322 are similar to pre viously 
monitored long-term cycles of solar-type star s ('see [Heni^ 
[19991 : [Lockwood et al.feoOTt [Hall et all [2009( 1. 

Figure 21 reveals that HD 130322’s brightness variabil¬ 
ity is anti-correlated with the strength of its H and K 
emission, as is common amo ng young, lower - main - 
sequence dwarfs. For example, [Lockwood et al.[ (|2007[1 
demonstrate the difference in behavior between young, 
solar-type stars with light curves dominated by dark 
spots and old solar-type stars with light curves dom¬ 
inated by bright faculae. In the young stars, photo¬ 


metric variability exhibits an inverse (or negative) cor¬ 
relation with chromospheric activity. In older stars, 
brightness variability and chromospheric activity exhibit 
a direct (or positive) correlation. Our Sun exhibits 
clear direct correlation between total sol ar irradiance and 
Ca H H and K emissio n (see Fig. 2 in [Lockwood et al.l 
[200711 . [Lockwood et al.[ ([2007[1 estimate the dividing line 



Fig. 4.— Top: The fluctuations in the Mt. Wilson S-index dur- 
ing the 13 observing seasons. Upper-Middle: Brightness of the 
primary target with respect to the Cl comparison star, measured 
with Keck I and the Til APT. Lower-Middle: Brightness of the 
primary target with respect to the C2 comparison star. Bottom: 
Differential magnitudes of the comparison stars, which show sta¬ 
bility to 0.0005 mag. The small variability in the two comparison 
stars, while the primary shows significant fluctuation, and perfect 
anti-correlation with seasonal mean brightness indicates that the 
variability in the light curves is intrinsic to HD 130322. The inverse 
correlation of S-index and brightness is typical of young, solar-type 
stars. 


ness variations to be around logi?^^ = —4.7. The 
origin al discovery paper of HD 130322b b y [IJdrv et al. | 
(200C) quoted a log R'^k value of -4.39 from [Santos et al.l 
([200II1 . corresponding to an age of only 0.35 Gyr and 
a rotation period of a pproximately 9 day s, according 
to the calibrations in (|Wright et al.l [20041 . This pre¬ 
dicted rotation period is much shorter than our observed 
Prnt = 26.53 days. By way of comparison. [Wright et al.l 
(j2004l l demonstrate that the correlation between chromo¬ 
spheric activity and stellar brightness in the G8 dwarf 
HD 154345 is positive, despite its similar properties to 
HD 130322. However, the logi?^^ value of HD 154345 
is -4.91 compared to -4.78 for HD 130322, which shows 
it to be much older than HD 130322 so that a positive 
correlation is expected. 


TABLE 4 

Summary of Photometric Observation eor HD 130322 


Observing 

Season 

(1) 

-^o6s 

(2) 

Julian Date Range 
(HJD - 2,400,000) 
(3) 

Sigma 

(mag) 

(4) 

Prot 

(days) 

(5) 

Full Amplitude 
(mag) 

(6) 

< P - Cl > 

(mag) 

(7) 

<P-C2> 

(mag) 

(8) 

< C2 - Cl > 
(mag) 

(9) 

2001 

99 

51912-52076 

0.0034 

23.0 ±0.2 

0.0046 ± 0.0009 

2.1189 ±0.0004 

0.5608 ± 0.0004 

1.5580 ±0.0001 
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TABLE 4 — Continued 


Observing 

Season 

(1) 

^obs 

(2) 

Julian Date Range 
(HJD - 2,400,000) 
(3) 

Sigma 

(mag) 

(4) 

Prot 

(days) 

(5) 

Full Amplitude 
(mag) 

(6) 

< P - Cl > 

(mag) 

(7) 

<P-C2> 

(mag) 

(8) 

< C2 - Cl > 
(mag) 

(9) 

2002 

230 

52288-52462 

0.0019 

24.0 ±0.2 

0.0022 ± 0.0003 

2.1193 ±0.0002 

0.5619 ± 0.0002 

1.5574 ± 0.0001 

2003 

79 

52645-52816 

0.0039 

29.0 ±0.4 

0.0067 ±0.0010 

2.1244 ±0.0005 

0.5664 ± 0.0005 

1.5580 ± 0.0002 

2004 

84 

53010-53189 

0.0042 

31.8 ±0.2^ 

0.0049 ±0.0012 

2.1253 ±0.0005 

0.5675 ± 0.0005 

1.5578 ±0.0002 

2005 

69 

53379-53551 

0.0029 

29.0 ±0.3 

0.0052 ± 0.0009 

2.1243 ±0.0004 

0.5668 ± 0.0004 

1.5575 ± 0.0002 

2006 

68 

53742-53913 

0.0020 

24.2 ±0.3 

0.0028 ± 0.0006 

2.1205 ±0.0003 

0.5625 ± 0.0003 

1.5580 ± 0.0002 

2007 

90 

54104-54282 

0.0025 

24.3 ±0.2 

0.0037 ±0.0006 

2.1195 ±0.0003 

0.5621 ± 0.0003 

1.5574 ±0.0002 

2008 

98 

54475-54637 

0.0034 

28.5 ±0.4 

0.0081 ±0.0007 

2.1254 ±0.0004 

0.5681 ± 0.0004 

1.5574 ±0.0001 

2009 

81 

54839-55003 

0.0042 

30.7 ±0.4 

0.0060 ±0.0012 

2.1261 ±0.0005 

0.5679 ± 0.0005 

1.5582 ±0.0002 

2010 

93 

55201-55382 

0.0069 

26.1 ±0.3 

0.0169 ±0.0014 

2.1278 ±0.0007 

0.5692 ± 0.0007 

1.5586 ±0.0002 

2011 

89 

55570-55738 

0.0031 

24.8 ±0.1=" 

0.0040 ±0.0008 

2.1229 ±0.0003 

0.5650 ± 0.0003 

1.5580 ±0.0002 

2012 

86 

55930-56098 

0.0019 

25.4 ±0.2 

0.0019 ±0.0006 

2.1200 ±0.0002 

0.5634 ± 0.0003 

1.5566 ±0.0002 

2013 

91 

56302-56470 

0.0046 

26.5 ±0.3 

0.0110 ±0.0008 

2.1249 ±0.0006 

0.5677 ±0.0005 

1.5573 ±0.0002 

2014 

99 

56659-56834 

0.0029 

27.7 ±0.2 

0.0042 ±0.0008 

2.1256 ±0.0003 

0.5681 ± 0.0003 

1.5576 ±0.0002 


^ Periodogram analysis gave half of the quoted period, sug¬ 
gesting that the star had spots on both hemispheres at those 
epochs. We doubled the photometric periods and their errors 
in these cases to get Prot- 

5. CONCLUSIONS 

Accurately determining the properties of planetary 
systems is extremely important as we move towards char¬ 
acterizing the thousands of exoplanets that are now con¬ 
firmed. It is only through understanding the host star 
that we are able to precisely measure the properties of the 
orbiting planet (s), which fuels both dynamic formation 
and evolution models. Through the Transit Ephemeris 
Refinement and Monitoring Survey (TERMS), we have 
studied HD 130322 because of the extensive RV coverage 
offered by HIRES, HRS, and CORALIE over the last 
^^14 years. The new and combined data has allowed us 
to determine a highly precise stellar radius of 0.85±0.04 
i? 0 , resulting in an updated Keplerian orbital solution to 
significantly limit the orbital dynamics of the b-planet. 
Through an extensive ^14 year photometric baseline at 
the APT, we carefully monitored the planetary phase 
during the predicted transit window, which did not re¬ 
veal any long-term variability of the host-star due to the 
presence of a companion. The HD 130322b planet had a 
transit probability of 4.7% at a depth of 1.57%. Signif¬ 
icant observations during the predicted transit window 
yielded a dispositive null result excluding a full transit 
to a depth of 0.017 mag and grazing transit to '^O.OOl 
mag. We were able to quantify the stellar rotation pe¬ 
riod with an unprecedented accuracy (26.53±0.70 days) 
by using the extensive photometric coverage. The vari¬ 
ation in differential magnitudes between the target and 


reference stars, as compared to the Mt. Wilson S-indices, 
also allowed us to better understand the stellar magnetic 
activity. However, the characterization of the HD 130322 
planetary system was only possible through the coming 
together of both collaborators and techniques, such that 
we were able to greatly improve the ephemeris of this 
system. The TERMS project consistently and systemati¬ 
cally provides accurate characterization of bright, nearby 
planetary systems, forwarding the understanding of exo¬ 
planets and their host stars in general. 
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